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ABSTRACT 


This report establishes a number of mathematical results 
relevant to the problem of constructing a triangulation, i.e., a 
simplicial tessellation, of the convex hull of an arbitrary 
finite set of points in n-space. 

The principal results of the present report are 

(a) A set of n+2 points in n-space may be triangulated in 
at most 2 different ways. 

(b) The "sphere test" defined in this report selects 
a preferred one of these two triangulations. 

(c) A set of parameters is defined that permits the 
characterization and enumeration of all sets of n+2 
points in n-space that are significantly different from 
the point of view of their possible triangularizations . 

(d) The local sphere test induces a global sphere 
test property for a triangulation. 

(e) A triangulation satisfying the global sphere 
property is dual to the n -dimensional Dirichlet 
tessellation, i.e., it is a Delaunay triangulation. 
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1 . Introduction 


Let P denote a set of m distinct points in n-space (n > 2, 
m > n + l.) Let C denote the convex hull of P . Assuming P does 
not lie entirely in some (n-l ) -dimensional manifold, we are 
interested in the problem of constructing a simplicial 
tessellation, T , of C such that each simplex, ter, has n+l 
points of P as its vertices, has non-null n - dimensional volume, 
and contains no other points of P . Such a tessellation will be 
called a triangulation. 

Another type of tessellation with which we shall be 
concerned is the Dirichlet tessellation. This is the set of. m 
cells, d , defined as 

d i = {q : \\q-p.i\ s Iiq-pyi for all Pj €P ) 

Thus the interior points of the cell d { are the points of n-space 
that are closer to the point p { than to any other point of p. A 
triangulation of P can be defined that is, in a certain sense, 
dual to the Dirichlet tessellation. Such a triangulation is 
called a Delaunay triangulation and will be defined in Section 9. 

Triangulations of point sets in 2-space are used in a 
variety of applications, particularly as an initial step in the 
analysis of data that is available at scattered points in the 
plane. After triangulation one may carry out interpolation, 
regridding to a rectangular grid, contour plotting, or other 
processes. The Dirichlet tessellation is used as a model in 
various scientific fields where it is appropriate to associate a 
unique region of space with each point in a finite point set. 
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Examples of applications of the Dirichlet tessellation are 
summarized and referenced in [Green '78]. 

Similar needs arise for the analysis of data defined at 
scattered points in higher dimensional spaces, but much less 
study has been devoted to triangulation algorithms and 
interpolation methods for these problems. 

Some early algorithms proposed for triangulation in the 
2 

plane required 0(m ) time. Since 1972 a number of computer 

programs for this problem, or the closely related Dirichlet 

tessellation problem, have been reported with estimates of 

13 15 

execution time in the range of m to m ' . See [Lawson '72], 
[Lawson '77], [Akima '78], [Green '78], [Cline ’84], and [Renka 
'84b]. Triangulation on the surface of a sphere is treated in 
[Lawson ’84] and [Renka ’84a]. The methods of various of these 
papers have been used in a number of proprietary graphics 
packages, and in the portable Fortran programs described in 
[Renka ! 84a] and [Renka ’84b] that are available from the ACM 
software distribution service. 

Bowyer [Bowyer '81] devised and implemented an algorithm for 
constructing the n-dimensional Dirichlet tessellation and the 
dual (Delaunay) triangulation. A very satisfactory execution 
time estimate of 0(a fc n ( ) was reported. The algorithm 
was implemented in ISO Fortran and performed well in a variety of 
test cases. 

Barnhill and Little [Barnhill ’84] presented ideas on a 
different approach to the n-dimensional triangulation problem and 
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also gave interpolation methods for use with triangular grids in 
3- and 4-dimensional space. 

The present report establishes some mathematical properties 
of n-dimensional triangulations that provide additional 
understanding of the problem. These results are generalizations 
of properties of the 2-dimensional problem given in [Lawson '77]. 

The algorithm given in [Lawson '77] operates by successively 
triangulating various 4-point subsets of the given point set P . 

It is easily seen that a set of 4 points in the plane admits of 
at most 2 different triangulations. In [Lawson '77] a local 
"circle test" was introduced that selected one qf these two 
possible triangulations. It was shown that the triangulation 
produced according to this criterion satisfied a global circle 
test property, a max-min angle property, and was dual to the 
Dirichlet tessellation. The equivalence 6f these latter two 
properties was established independently by Sibson in [Sibson 
■78]. 

The principal results of the present report are 

(a) A set of n + 2 points in n-space may be triangulated in 
at most 2 different ways. 

(b) The "sphere test" defined in this report selects 
a preferred one of these two triangulations. 

(c) A set of parameters is defined that permits the 
characterization and enumeration of all sets of n+2 
points in n-space that are significantly different from 
the point of view of their possible tr iangularizations . 

(d) The local sphere test induces a global sphere 
test property for a triangulation. 
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(e) A triangulation satisfying the global sphere 

property is dual to the n- dimensional Dirichlet 
tessellation, i.e., it is a Delaunay triangulation. 

2 . Barycentric coordinates and their geometric interpretation 
The convex hull of n+1 points is called a simplex. The 
convex hull of a subset consisting of n of these points is called 
a facet of the simplex. 

Let p , ..., p n be n+1 distinct points in n-space. 

Define the matrix 

s = [ 1 1 * * ’ 1 ] 

L P 1 p 2 • • • P n + 1 J 


Let t be the (possibly degenerate) simplex with vertices p^, 
. .., . The n - dimensional volume of t is given by 


Vol(t) = | Det ( B ) | /n ! 

Let q be an arbitrary point in n-space. If t is 
nondegenerate, i.e., if Vol(t) * 0, the numbers, b , . . . , , 

satisfying 


(l) 


L p l P 2 


1 I 


n + 1 J 


1 b 1 

L «+l J 



1 

I 

J 


are called barycentric coordinates of q relative to the simplex, 
t . 
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For each s, the sign of £> s indicates the position of q 

relative to the hyperplane, , containing the facet of t 

opposite vertex p . Thus b =0 when q is in H , b > 0 when q 

s s s s 

is on the same side of s as p , and b < 0 when q is on the 

3 3 S 

opposite side of from p g . 

Some consequences of these facts are: 

(a) The point, q, is in the simplex, t, if and only if all 

b >. 0. 
s 

(b) If q is strictly outside the simplex, t, then one or 
more of the *> s 's are negative. These negative b s 's 
identify the facets of t with whose vertices q can be 
connected to form nondegenerate simplices neighboring 
to t. 

(c) At least one of the fc g 's must be positive since 

Zb = 1 . 
s 

3 . Trianqulations of sets of n+2 points 

In the 2D problem, study of all possible ways of 
triangulating the convex hull of four points provided the key 
ideas that led to a systematic procedure for improving the 
triangulation of the convex hull of larger sets of points as one 
point at a time was introduced. The study of sets of n+2 points 
provides similar insights in the n-dimensional problem. 

In 2-space all of the different possibilities for 
configurations of four distinct points, not on a common line, are 
illustrated by the three cases in Figure 1. 


5 



CASE 1 




Figure 1. Configurations of four distinct 
non-colinear points in 2-space. 


In Case 1 there is a unique triangulation using two 
triangles. In Case 2 there are two possible triangulations, each 
using two triangles. In Case 3 there is a unique triangulation 
using three triangles. 


One might expect the general n-dimensional case to be much 
more complicated than the 2-D case. Although there is more 
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complexity in the higher dimensional cases, there are, 
fortunately, two very useful properties that persist. It remains 
true that a given set of n+2 points in n-space admits of at most 
two different triangulations by n-dimensional simplices. 
Furthermore, if the n+2 points do not all lie on a common 
n-sphere, one of the possible triangulations is uniquely selected 
by the sphere test that will be described shortly. 

These facts will be established in the following two 
sections. 


4. Signature sets 






Theorem 1 

Let 

P be a 

set of n+2 

point Si p , . 

... P n+2 . in 

n-space not. lying entirely 

In any (n- 

l)-di mensi ona I 

mani fold . 

There 

Is a partitioning of 

the index 

set {1, 2. 

. n+2 } into 

three 

sets. Sq , 

V 

and S 2< 

and a set 

of numbers, c 

. . sat i sfy i ng 

(2) 



Z C D 

/€S 1 

» X 

c / p / 


( 3 ) 



X C . 

'«! ' 

“ * C / 
,€S 2 

= 1 


(4) 



c ■ . 0 , 

' 6 s o 



(5) 



c . > 0. 

t 

i € S,US rt 
1 2 



The numbers c. are uniquely determined by 

sets S^. s „ , and S_, are also unique, with the 
0 12 ^ 

the labeling of and S ^ could be arbitrarily 
Proof. The (n+l)x(n+2) matrix 


the set P . The 
understanding that 
interchanged . 


A = 



1 

P n+2 


1 

I 

J 
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is of rank n + 1 due to the hypothesis that the points of P do not 
all lie in any (n-l ) -dimensional manifold. The system 
( 6 ) Ax - 0 

* 

has a 1-dimensional space of solution vectors. Let x be a 

nonzero vector satisfying Eq.{6) and normalized to satisfy 

2|x. | = 2. Since the first row of Eq.(6) is zx - 0, the vector 
* 

x must have both positive and negative components, and the sum 

of the positive components must equal the sum of the magnitudes 

of the negative components. In fact, with the specified 

normalization, these sums must each be 1. 

* 

If any components of x are zero, let s Q be the index set 

for these components. Let 5 1 be the index set for components of 

* * 
x of one sign and let be the index set for components of x 

of the other sign. Both and s 2 are non-null. Define 

♦ 

o , I x . } f i — 1# ••• / n+2 

Then Z ( c . = 2 and all of the Eqs.(2-5) are satisfied. 

This specification of the sets, s Q , and s , and the 

numbers, c . , is unique to within the possible interchange of Sj 

and S 2 , since if Eqs.(2-5) were satisfied by any other sets, S' Q , 

s ^ , and s * , and numbers, c‘. , it would permit the construction of 

a solution vector for Eq.(6) not in the same 1-dimensional 
* 

subspace as x . □ 

For any point set P satisfying the hypotheses of Theorem 1 , 
the associated sets, S Q , s , and s 2 , will be called the signature 
sets of p . 

* 

Let p be the point given by the equal left and right 

* 

members of Eq.(2). Note that p is in the convex hull of 
{ p.ties and also in the convex hull of {p j :deS 2 >. In fact the 
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geometric interpretation of Theorem 1 is that the subsets of P 

indexed by and s 2 are. the smallest two disjoint subsets of P 

whose convex hulls have a point in common, and the common point, 

* 

p , is unique . 

As examples, in Fig. 1, we can identify the signature sets 

* 

and the common point p as follows: 

Case I S Q = {2}, S x = (4), S 2 = {1, 3), p* = p 4 

Case II s q = Null , S = (1, 3}, S = { 2 , 4}, p = the 

intersection of lines P^P 3 and. p p 4 « 

Case. Ill S Q -= Null, S 1 = {4}, S 2 - {1, 2, 3>, p* = p 4 

To relate this to the possible triangulations of the convex 
hull of P, we next identify the possible nondegenerate simplices 
that can be formed from the points of P . 

For each / = 1, ..., n+2, let denote the (possibly 
degenerate) simplex formed using the points P\p as vertices. 

(The notation /’Np^. denotes the subset of p remaining when point 
p { is removed.) The n- dimensional volume of is Ac { where the 
's are given by Theorem 1, and A is a positive constant 
independent of / . Thus the nondegenerate simplices are just the 
set (t. : ies^us^. 

What are the possible groupings of these simplices to form a 
triangulation of the convex hull of P? 

5 . Admissible triangulations of n+2 points 

Theorem 2 Let P be a point set as in Theorem 1. There are 

at most two distinct tr i angu l at i ons of the convex hull of P. 

namely, T, - (t.r/es,) and T n where the sets S, and 

3 1 i 1 2 i 2 1 
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S 2 are as defined in Theorem 1, and t . is the simplex with vertex 
set P-{p.}. If one of the sets or S 2 is of cardinality 1-. the 
correspond i ng set T or T 2 is not an adml ss i bl e triangulation. 

Proof . This theorem will be proved by establishing the 
following four assertions: 

(a) Any pair of simplices, one indexed in and the 
other indexed in s^, is mutually overlapping and 
thus cannot be used in the same triangulation. 

(b) Any pair of simplices, both indexed in or both 

indexed in is nonoverlapping, and thus the 

sets r 1 = {t.-.ias^} and T 2 { t . : i es > are each 
nonoverlapping sets of simplices. 

(c) Each of the sets r and T ^ covers the entire 
convex hull, c, of f, i.e., any point q ec is also 

«r 

contained in some and in some t^er 2> 

(d) T. is not an admissible triangulation of c if 


|s. | = l. 

Note that Eq.(2) can be solved for any one of the p. ' a , 
i es ^us 2 , and the resulting equation gives the barycentric 
coordinates of one point of P with respect to the simplex formed 
by the others. For example, choose an index, kes^, and solve 
Eq.(2) for p fc , obtaining 

(7) p R = z <c /c Ip - I (c ( /c )p + £ (c /c Ip 
/es 2 iaS Q 

from which we may write the barycentric coordinates of p 

A 

relative to the simplex t as 

A 


(S) 

(9) 

( 10 ) 


b. = c./c, > 0 for i as n 

i i k 2 

b. = -c./c, <0 for /«(S -{*}) 

i i k ' 1 

b. =0 for i as 

i 0 

- 10 - 



Eq.(8) implies that for each /eSg, p. and p k are on the same 
side of the common facet shared by t and t . Thus t. and t 

t fC / ft 

overlap, establishing assertion (a). Eq.(9) implies that for 

each ies^Xk, p. and p k are on opposite sides of the common facet 
shared by t and Thus t. and t k do not overlap, establishing 

assertion (b). 

Let q be an arbitrary point in c. Then there are 
non-negative coefficients, d . , such that 

n+2 

(11) q = I dp 

i=l 

n+2 

(12) 1 = Z d . 

/ = 1 ' 

Rewrite Eqs.(2-3) as 


(13) 0 = 2 c p - 2 c p. 

, » i i / / 

/es x /€S 2 

(14) 0 = 2 c - 2 c . 

/€S 1 /€S 2 1 

Using an indeterminate, A, form A times Eq.(13) plus 
Eq.(ll), and A times Eq.(14) plus Eq.(12): 


(15) 

(16) 


2 (d.+Ac.)p. + 2 (cf.-Ac.)p. + 2 d . p . 


2 (cf.+Ac.) + 2 d.-Ac.) + 2 d. 

. „ / / . ' / / . _ / 
/6S, /€S. / € S _ 

1 2 0 


There is a range of values of A , A < A < A , for which 

man max 

all of the coefficients, (d .+ Ac. ) , ,/es , and (d.-kc.), ies^, 
appearing in Eqs. (15-16) are nonnegative. At the low end of this 


range, at least one of the coefficients indexed in say 

dj+hCj, is zero, showing that qety Similarly at the high end of 
the A range at least one of the coefficients indexed in s , say 
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cf - Ac^, is zero, showing that q These limiting values of A 

are given by 


A . = max { -c 1 . /c . : i e s } 
min 'ii 1 

A m = v = m±n{cf /c :/es } 
max / / z 

This establishes assertion (c). 


Consider now the case in which one of the sets S 1 or s 2 has 

cardinality 1 , e.g. , suppose = 1. Then the triangulation, 

T , consists of a single simplex, say t . . The point p. that is 

J J 

not a vertex of t . is therefore not a vertex in the 


J 

triangulation, r . Thus the triangulation r ^ is not admissible 
because it does not include all points of P as vertices. □ 


In a different context, namely in developing a stable method 
for evaluation of multivariate splines, Grandine [Grandine '84] 
proved a theorem encompassing a subset of the above Theorem 2 . 

His theorem states that an arbitrary point in the convex hull of 
a set of n +2 points in n-space can be in the interior of at most 
two of the n +2 simplices that can be formed using these points. 

As examples of Theorem 2, consider again Figure 1. In Case 
1, with s = {4} and S = {1, 3}, the triangulation T = { £ ,» 

X 4 Z i. o 

is admissible while = ( £ 4 ) is not, again because p ^ is not 

used as a vertex. In Case 2, with = { 1 , 3} and S = {2, 4}, 

there is a choice of two admissible triangulations, T = {t , t } 

2 ± si 

and r 2 = {t 2 , £ 4 >- In Case 3, with S 1 = {4} and S = {1, 2, 3), 
the triangulation T = ( £ lf £ 2 ' 1 3 ^ is ad ®issible while T = { £ 4 > 
is not, because p is not used as a vertex. 

Corollary 1 In the context of building tr i anqu 1 at i ons , an 
enumeration of all possible s i gn i f i cant l y different 
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configurations of n+2 distinct points in n- space . not lying in 
any C n-li -d i mens i onal manifold, is given, by all of the possible 
ways of assigning values to | S Q |» |S^|, and | S g ( satisfying 

(17) jS 2 | 2: |S 1 | 2: 1 

(18) • |S 2 | > 2 

(19) |S Q | > 0 

(20) |S Q | + I S ± ( + |S 2 | = n + 2 

Proof. The sets and S 2 must each be nonempty to satisfy 
Eq.(3). Since the sets S 1 and s 2 are not mutually distinguished 
we may arbitrarily use s 1 as the label of the smaller of the two 
sets when they are not of the same size. These considerations 
establish Eq.(17). 

The sets and S cannot both be singletons for then they 
would have to index the same single point in order for Eq.(2) to 
hold. This is ruled out by the hypothesis that all points of p 
are distinct. Along with the convention of Eq.(17) this gives 
Eq (18). 

Eqs. (19-20) follow from previous discussions. Q 

Using Corollary 1 one may list the values of (l s 0 l» |'S | , 

| s 2 | ) for all significantly different configurations of n+2 
points in n-space. This is done for dimensions 1, 2, 3, and 4 in 
Table 1. 
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Table 1. Characterizing parameters for all of 
the significantly different configurations 
of n+2 points in n-space for n = 1, 2, 3, and 4. 

n=l n = 2 n = 3 n = 4 


I S 1 I i S 2 ' 


I s n I l s J l S ol I s n I 


1 

0 

0 


1 

2 

1 


2 

2 

3 


2 

1 

1 

0 

0 


' S ll 


\S. 


1 

2 

1 

2 

1 


2 

2 

3 

3 

4 


3 

2 

2 

1 

1 

0 

0 

0 


1 

2 

1 

2 

1 

3 

2 

1 


|s 2 l 


2 

2 

3 

3 

4 

3 

4 

5 


For each dimension in Table 1 the triples are listed in 
reverse lexicographic order. With this ordering cases in the 
same row have the same pair of values of | s ^ | and | s 2 | , differing 
only in |S Q | . Since the possible triangulations are determined 
by | s ^ j and | s 2 | , there is a significant geometric relationship 
between the possible triangulations for cases appearing in the 
same row. This will be explained further after introducing 
Figure 2. 

The three cases shown for n - 2 are those previously 
illustrated in Figure 1, and in the same order. Diagrams 
illustrating the five cases for n - 3 are given in Figure 2. 
Recall that cases with |S 1 | = 1 admit only one triangulation each 
while cases with | s ^ | > 1 admit two distinct triangulations each. 
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non-coplanar 


Note that in the first three cases of Figure 2 the base 

facet of each diagram has the same configuration as the 

corresponding case in Figure 1. Case 1 admits of a further stage 

of reduction since the base edge in Figure 1 is the diagram for 

the case of 3 points in 1-space, i.e., the case corresponding to 

the first and only entry for n = 1 in Table 2. 

These are examples of the following general principle: If a 

point set P has parameters (o Q a^, a ) with o Q * 0, the subset 
* 

p consisting of the points indexed in and will lie in an 

(n-o Q ) -dimensional manifold and have a (0, o^, a configuration. 

(Note that these two related configurations are in the same row 

in Table 2.) Any triangulation of P is related to a tesselation 
* 

of p by the fact that any simplex in a triangulation of P is the 

* 

convex hull of the union of a simplex in a triangulation of P 

* 

with the points of p\p . 

6 . Configurations having all pairs of points connected 

Note that in Case 3 and in the second triangulation shown 
for Case 4 in Figure 2 , every pair of points is connected by an 
edge used in the triangulation. This is true of any 
triangulation of n+2 points in n -space if the triangulation 
consists of three or more simplices. This follows from the 
observation that if any three different subsets of size n + 1 are 
selected from a set of n+2 points, every possible pairing of 
points must appear in one or more of the selected subsets. 

This observation has. particular significance in dimensions 4 
and higher, since then there are configurations of n+2 points 
having 3 or more points in both and S ; e.g. , see Row 6 for 
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n = 4 in Table 1. In such cases both of the possible 
triangulations will have all points connected by edges of the 
triangulations. Therefore the two possible simplicial 
triangulations are not distinguished from each other by 
information about connectivity of pairs of points. 

We shall return to this point in Section 9 in connection 
with the Dirichlet tessellation. 

7 . The sphere test for a set of n ± 2 points in n-soace 

Let P continue to denote a set of n + 2 points in n-space, not 
lying entirely in any n-1 dimensional manifold. In the preceding 
section it was seen that there may be either one or two ways to 
triangulate the convex hull C of such a set P . For the cases in 
which two triangulations are possible we introduce the sphere 
test as a way of choosing one of the triangulations. 

Suppose r is a configuration that admits two possible 
triangulations. Using the notation and results of Theorems. 1 and 
2, it follows that |S 1 | > 2, | S 2 1 i 2, and n i. 2. The two 
possible triangulations are T ^ = {t . :ies } and r 2 = .{t { :ies 2 >, 
where t { is the nondegenerate simplex with vertex set / > \p / . 

For each ies^JS^, let E j be the unique n-sphere 
circumscribing the simplex and let be the open n-ball whose 
boundary is £. . 
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The Sphere Test 


If all points of P lie on the same sphere, i.e., 

all of the spheres E. are coincident, then the 

sphere test does not distinguish between and r^. 

Otherwise choose /es,us_. If p ,«B . select the 

12 r jj 

triangulation T or T that includes t , while if 

* ^ y 

p ,€S . select the triangulation T or T that does not 
J J 12. 

include t . . 

J 


As an example consider Case 2 of Figure 1 with = {1,3} 
and S g = { 2 , 4}. To apply the sphere test we may choose any one 
of the four points and ask whether it is inside or outside the 
open ball circumscribing the other three points. Figure 3 
illustrates the four possible ways of applying the test for this 
example. We find p eB , p.es p *B , and p es r Thus any one 
of these tests results in the selection of r 2 = {* 2 ' * 4 ) as the 
preferred triangulation. 
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Figure 3. Example of application of the sphere test. 

The crucial fact that the result of the sphere test is 
unique independent of the choice of the test point is established 
by the following theorem. 

Theorem 3 Let P be a point set satisfying the hypotheses of 
Theorem 1 with signature sets satisfying | S [ £ 2 ancf l s 2 l - 2 . 

If the points of P do not all lie on the same n-sphere then 
either 


P/**/ 

for 

al i 

“*1 

and 

p eB . 

r i i 

for 

al l 

/eS 2 

p/“/ 

for 

al l 

lts 2 

and 

p.eB 

r i i 

for 

al l 

,€5 1 
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For the proof of this theorem it is useful to have the 

following two lemmas on intersecting spheres. 

Lemma 1 Let E and E be two distinct n-spheres in n-space, 

(n > 2), intersecting in a set E ^ consisting of more than a 

single point. '■ Then 13 an (n-D-sphere contained in a 

uniguel y determined hyperplane 

Lemma 2 Let B^ and B 2 denote the open balls bounded by E ^ 

and E respect i ve 1 y . On one side of the hyperplane, B^ 

will contain B , while the reverse inclusion will prevail on the 

other side of h^ 2> Let denote the open halfspace on the side 

of ft , _ in which B . contains and let denote the other open 
12 1 2 2 r 

halfspace. Then 

H OB c H ,r>B , 

12 11 

(ff 2 nB 2 ) n (« 1 n£ 1 ) = Null 
b 2 n (« 1 n£ 1 ) = Null 

(24) H 2 R B i c « 2 nB 2 

(25) ( H 2 nS i ) n ( H 2 n£ 2 ) " Nul1 

(26) B 1 n ( H 2 n B 2 ) = Null 

Proof of Lemmas . The validity of these lemmas for n=2 is 
clear from consideration of Figure 4. 


( 21 ) 
( 22 ) 
(23) 
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Figure 4. Configuration of Lemmas 1 and 2 in 2-space. 

In higher dimensional spaces the configuration of the two 
intersecting balls is symmetric about the line, I, connecting the 
centers of the two balls. Thus the intersection of the objects 
E i # E g , Sg, h 12 , H and « 2 with any 2-dimensional manifold 

containing t again gives the configuration of Figure 4. We omit 
further details of the proof. □ 


The principal conclusions to be used subsequently from Lemma 
2 are Eqs.(23) and (26). For example from Eq.(23) we know that 
if a point, q, is in then q«£ 2 * 


Proof of Theorem 3. Without loss of generality let jes and 
suppose p jtB . so that T 1 is selected. It suffices to show that 
p .€B . for every and p .eB . for every ies^. 

Let ies^JS^. Let t denote the (n-1 ) -simplex forming the 
common facet of the simplices t and t . Thus t is the 

J J 

(n-l ) -simplex with vertex set P\{p.,p.), 
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Let h.j denote the unique hyperplane containing t . The 
intersection of the n-spheres E. and Ej is the ( n -1) -sphere E 


that circumscribes t 


... E . . lies in ft . . 
ij ‘J U 


On one side of ft . . , B . 

‘J i 


is contained in B . , while on the other side of ft . . , B . is 

J U J 

contained in B . . From information about one point of E, , namely 

Pj , that is not in ft the relative containment relations 

between B. and B . can be determined. 

• J 

Let H j denote the open half space on the same side of ft as 
Pj and let W J denote the open halfspace on the other side of ft_ . 
See Figure 5 . 



Figure 5. Configuration of Theorem 3 in 2-space. 


Since p ,<B . but p .e£ . , it follows that 
J J J ! 

( 27 ) B r\H . C b n« , 

J J I J 

Then, as in Lemma 1, in the other half space, H j , we have 

( 28 ) a.nHl c 

‘ J J J 

We now consider the two possible cases of or 

If , € S 1 then p i and Pj are on opposite sides of ti . and thus 


22 


p.eH'j. Also p .eE j and thus p jf\H ' . From Eq.(28) and Lemma 2 
it follows that p . . 

If i€S 2 then p { and Pj are on the same side of ft and thus 
p^Hj. Also Pj* E j and thus p^EjCUij. From Eq.(27) it follows 
that p.eB.. □ 

8 . The local and global sphere tests 

Let T be a triangulation of a point set p in n-space. The 
triangulation T satisfies the global sphere test if for each 
simplex t the open n-ball B. circumscribing t contains no 
points of P . 

A pair of n-simplices, t { and t , sharing a common 
(n-1 ) -dimensional facet t^j will be said to satisfy the local 
sphere test if the vertex of t. not in t is outside the open 
ball B j circumscribing t y By Theorem 3, this is equivalent to 
the requirement that the vertex of t not in is outside the 

open ball s / circumscribing t.. 

Theorem 4 If a triangulation T of a point set P has the 
property that every pair of simplices shar t ng a common facet 
satisfies the local sphere test, then T satisfies the global 
sphere test . 

Proof . The proof will be by contradiction. Suppose the 

* 

hypothesis is satisfied, but there is some point p ep and some 
simplex t'er such that the open ball B' circumscribing t* 

* 

contains p . 

• * 

Let i be a line segment from p to some point q interior to 
t‘ . By a small perturbation of the position of the end point q, 
if necessary, we may assume that wherever i passes from one 
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simplex to another it intersects the relative interior of a 
common facet. 

Relabel the simplices intersected by i so they are denoted 

by t , t 2 , , ordered along the line i from which was 

* 

previously called t' , to t , which has p as one of its vertices. 

For i = 2, . .., k , let p { . be the vertex of t. that is not a 

vertex of t . . For / = 1, ...» k , let E. be the sphere 

circumscribing and let B. be the open ball whose boundary is 

* 


E . . Note that p, = p and B , = B ■ . 

/ k 1 

By hypothesis, p ,£B . _ , / =2, ..., K, but we are assuming 
p ea . Since p eB and p <B , there exists a smallest index, 

X X n 1 K K X 

j , such that p.€B . and p «fl . Let h denote the unique 

k j k j + i 

hyperplane containing the facet common to t . and t + ^ . Figure 6 
illustrates t^, t j + ^> and related objects for the case of n = 2. 



Figure 6. Configuration of Theorem 4 in 2-space. 


Let H ^ be the open halfspace on the same side of h as 
t j + l’ Let H j be the opposite open halfspace, i.e., the halfspace 


on the same side of h as t 


Then 


and 


p . , e £ n H . , 
y J + 1 7 + 1 7 + 1 


p . , e b . 
7+1 7 
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thus 


Wi c Vi n V i 


p„ € B ,C\H . , , 
** 7 7+1 


p, e S . ,0H . „ 

** 7+1 7+1 


contradicting the assumption that p^.«S^ + 1 . 0 


9 . The global sphere test and the Dirichlet tessellation 

Let P again denote a finite set of distinct points in 
n-space, not lying entirely in any (n-1 ) -dimensional manifold. 
With each point, p .ef , we associate the cell, 

a t = {q : llq-pyi S llq-pyi , for all Pj ep ) 

The cell d. is called the Dirichlet cell associated with p t 
relative to the point set P. ' The set of all d^'s is called the 
Dirichlet tessellation of n-space associated with the point set 
P. Clearly the cells d are disjoint except for common 
boundaries, and the union of all of the d^'s covers all of 
n-space. 

Each cell, d . , is the intersection of a finite number of 
half spaces. Each such half space is bounded by the hyperplane 
that perpendicularly bisects the line segment connecting the 
point p i to another point of P . 

Let Q denote the set of points, . , that occur as vertices 
of the d.'s. Each q . is the unique intersection point of at 
least n facets of some cell, d., and thus is the unique 
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intersection point of at least n of the bisection hyperplanes, 
with some n of them being linearly independent. 

Let r . = llq ,-p .ii . Since q .ecf . , there are no points of P 
whose distance from q . is less than r j . For each bisecting 
hyperplane on which g^. lies, there is another point, pep, which 
is the reflection of p. relative to this hyperplane and whose 
distance from q^ is also r .. Thus there are n+1 or more points 
of P at the distance from q^ , and this set of points does not 
lie in any (n-l ) -dimensional manifold. 

Conversely it can be verified that any point in n-space that 
attains its minimum distance from points of P at n+1 or more 
points of p that do not lie in an (n-l ) -dimensional manifold must 
be a vertex of one or more of the cr^ *s, i.e., must belong to the 
set Q . 

With each point q^€Q , associate the convex hull, , of the 
points of P that are at the minimal distance, , from q^ . The 
set, s ’ , is called a Delaunay cell. The set of all s^.'s 
constitutes the Delaunay tessellation of the convex hull of p . 

In particular it can be verified that the s^.'s are mutually 
disjoint except for common boundaries, and the union of all of 
the s^'s coincides with the convex hull of P. 

The Delaunay tessellation and the Dirichlet tessellation are 
dual to each other in the sense that each cell, d . , of the 
Dirichlet tessellation is associated with a vertex, p , of the 
Delaunay tessellation, and each cell, , of the Delaunay 
tessellation is associated with a vertex, q^ , of the Dirichlet 
tessellation. 
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If the points of P are in "general" position, each cell, s , 
will have just n+1 vertices, and thus will be a simplex. In this 

case, the Delaunay tessellation may be called a triangulation. 

In practical applications, e.g., [Bowyer ’81], where one uses the 
Delaunay tessellation as a means toward producing a 
triangulation,, one can replace any cell, s^, that has more than 
n+1 vertices by an arbitrary triangulation of that cell, thus 
producing an overall triangulation of P. It appears common to 
extend the name Delaunay tessellation to such a triangulation. 

We may now observe that a triangulation satisfying the 
global sphere test of Section 8 is, in fact, a Delaunay 
tessellation, possibly in the extended sense just mentioned. Let 
T be a triangulation of P satisfying the global sphere test. 

With each simplex, ter, associate the point, q, at the center of 

the circumsphere of t . If two of more simplices have the same 
circumcenter , replace these simplices by their union, s. Note 
that all vertices of such a cell, s, lie on a common sphere. 

The circumcenter points, q, associated with the cells of 
this tessellation satisfy the properties of the set, Q, noted 
previously. Thus these cells are all Delaunay cells. Since 
their union covers the convex hull of P, no Delaunay cells are 
missing, so this is a Delaunay tessellation for P . 


27 



References 


Akima, Hiroshi (1978), A method of bivariate interpolation and 
smooth surface fitting for irregularly distributed data 
points, ACM Trans. Math. Software, 4 , 148-159; also 
Algorithm 526, 160-164. 

Barnhill, R. E. and Little, P. F. (1984), Three- and 

four-dimensional surfaces, The Rocky Mountain Jour. Math., 
14, 77-102. 

Bowyer, A. (1981), Computing Dirichlet tessellations, The 
Computer Journal, 24 , 162-166. 

Cline, A. K. and Renka, R. L. (1984), A storage-efficient method 
for construction of a Thiessen triangulation, The Rocky 
Mountain Jour. Math., 14 , 119-140. 

Grandine, T. A. (1984), The stable evaluation of multivariate 

B-splines, MRC Technical Summary Report No. 2744, University 
of Wisconsin Mathematics Research Center, Madison, 

Wisconsin, 14 pp. 

Green, P. J. and Sibson, R. (1978), Computing Dirichlet 

tessellations in the plane. The Computer Journal , 21 , 
168-173. 

Lawson, C. L. (1972), Generation of a triangular grid with 

application to contour plotting, Jet Propulsion Laboratory 
Internal Technical Memorandum No. 299, Pasadena, California. 

Lawson, C. L. (1977), Software for C 1 surface interpolation, in: 
J. R. Rice, ed.. Mathematical Software III, Academic Press, 
New York, 161-194. 

Lawson, C. L. (1984), C 1 surface interpolation for scattered data 
on a sphere. The Rocky Mountain Jour. Math., 14 , 177-202. 

- 28 - 



Renka, R. J. (1984a), Interpolation of data on the surface of a 
sphere, ACM Trans. Math. Software , 10, 417-436; also 

Algorithm 623, 437-439. 

Renka, R. J. (1984b), Algorithm 624. Triangulation and 

interpolation of arbitrarily distributed points in the 
plane, ACM Trans. Math. Software , 10, 440-442. 

Sibson, R. (1978), Locally equiangular triangulations. The 
Computer Journal , 21, 243-245. 


29 



